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ABSTRACT 

Many classical nova ejecta are enriched in CNO and Ne. Rosner et al. recendy suggested that the enrichment 
might originate in the resonant interaction between large-scale shear flows in the accreted H/He envelope and 
gravity waves at the interface between the envelope and the underlying C/0 white dwarf. The shear flow 
amplifies the waves, which eventually form cusps and break. This wave breaking injects a spray of C/0 into 
the superincumbent H/He. Using two-dimensional simulations, we formulate a quantitative expression for the 
amount of C/O per unit area that can be entrained, at saturation, into the H/He. The fraction of the envelope that 
is enriched depends on the horizontal distribution of shear velocity and the density contrast between the C/O 
white dwarf and the H/He layer but is roughly independent of the vertical shape of the shear profile. Using this 
parameterization for the mixed mass, we then perform several one-dimensional Lagrangian calculations of an 
accreting white dwarf envelope and consider two scenarios: that the wave breaking and mixing is driven by the 
convective flows; and that the mixing occurs prior to the onset of convection. In the absence of enrichment prior 
to ignition, the base of the convective zone, as calculated from mixing-length theory with the Ledoux instability 
criterion, does not reach the C/O interface. As a result, there is no additional mixing, and the runaway is slow. 
In contrast, the formation of a mixed layer during the accretion of H/He, prior to ignition, causes a more 
violent runaway. The envelope can be enriched by < 25% of C/0 by mass (consistent with that observed in 
some ejecta) for shear velocities, over the surface, with Mach numbers < 0.4. 

Subject headings: hydrodynamics — nuclear reactions, nucleosynthesis, abundances — novae, cataclysmic 
variables — methods: numerical — waves 



1. INTRODUCTION 

Classical novae are a manifestation of thermonuclear run- 
aways in accreted hydrogen shells on the surfaces of white 
dwarf stars in clo se binary systems (see, e.g., the review by 
iGehrz etaOllQQSh . Compelling observational data indicate 
that the material ejected by some classical novae can be sig- 
nificantly e nriched i n C, N, O, and Ne, by > 30% by mass 
JCjehrzetaL.1998.: l.ivio & Truranltl994) . Since the abun- 
dance of CNO catalysts constrains the rate of energy release, 
such high levels of CNO enrichment are required in the fastest 
novae (Truran 1982), for which the hydrogen burning reac- 
tions (via the CNO cycle) increase the entropy of the enve- 
lope faster than it can adjust. It was recognized early that if 
some of the underlying white dwarf matter"' could be mixed 
into the accreted envelope prior to the final stages of the ther- 
monuclear runaway, then the explosion would be more en- 
ergetic and the ejecta more enriched in CNONe elements 
jStarrfield et alJIT972l) . One-dimensional models that best 
reproduce observations typically accrete material "seeded" 
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with a super-solar c omposition (e.g., iHernanz et alJ Il996t 
'Starrfield et alJl998^ . 

The question of how the accreted envelope is enriched 
has challe nged theory for sever al decades (see, e.g. the 
review by tLivio & Truranl[T99l . For very slow accretion 
(< lO^^^MQyr^'), the downward diffusion of H into the 
underlying C/O or O/Ne layers ("Kovetz & Pri alnikI 119851 
Prialnik & Kovetz 1984) could trigger ignition in the H 
diffusive tail and drive convective mixing of heavy elements 
into the envelope during the early stages of runaway. It is 
unclear, however, whether there is time at higher accretion 
rates for this process to be relevant. The possibility of mixing 
via convective overshoot was first considered by Woosley 
( 1986). Numerical simulations of convective penetration 
and mixing at the core-envelope interface have been car- 
ried o u t in both two dimensions ( Glasner. Livne. & TruranI 
199?; 'Kercek, Hillebrandt, & Truratt 11991 and three 
( Kercek, Hillebrandt, & Truran 1999fc. From comparisons 
of two- and three-dimensional simulations and a careful 
resolution study, Kerce k et al. ( .199,8^ 1999J concluded that 
convective penetration would not significantly enrich the 
accreted envelope. 

The general problem o f she a r mixing was consid- 
ered by m any ( Fuiimotg| 119881 iKippenhahn & Thomas I 
Il978t iMac Donald 1983,), but to date no self-consistent 
calculation s of these effect s have been performed. 
iKippenhahn & Thomas I (119781 ) considered shearing from 
differential rotation in the envelope and estimated the amount 
of mass mixed by assuming that the envelope was marginally 
stable according to the Richardson criterion. MacDonay 
i 19831) considered non-axisymmetric perturbations and the 
resulting redistribution of matter and angular momentum. 
Fuiimoto (1988) emphasized the role of barotropic and 
baroclinic instabilities in transporting angular momentum 
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through the envelope. A vari ation on this was envisaged by 
Uparks & Kutter ( 19|^ and iKutter & SpaiQ ^99,9), who 
suggested that convection just prior to the runaway would 
transport angular momentum inward and lead to a large 
horizontal shear above the C/0 white dwarf. 

While these studies examined the generation of shear in the 
envelope, the details of how C/O is mixed into the envelope 
and the quantity of mass mixed, as functions of the shear pro- 
file, were not resolved. In this paper, we examine a promising 
mechanism for effecting the entrainment: a resonant interac- 
tion between large-scale shear flows in the accreted envelope 
and interfacial gravity waves (^losner et al. 2001 ). The greater 
compositional buoyancy in the C/O white dwarf means that 
the interface sustains gravity waves. Miles ( 1957) showed that 
in the presence of a shear flow (i.e., a "wind"), gravity waves 
with a group velocity matching a velocity in the shear flow 
are resonantly amplified. These waves eventually form cusps 
and break. When the waves break, they inject, analogously to 
ocean waves, a spray of C/O into the H/He atmosphere. 

The source of the shear could arise from a number of mech- 
anisms, such as those described above. A self-consistent cal- 
culation of the accretion boundary layer and the wind profile 
in the atmosphere is far beyond the scope of current simu- 
lations. Indeed, as we explain in § |2] the resolution needed 
to simulate the interaction is so fine as to prohibit the direct 
simulation of this effect within a macroscopically large (i.e., 
the entire atmosphere of the white dwarf) computational do- 
main. What this paper provides is a scaling relation (eq. 13) 
for the mixing rate per unit area, for a given wind velocity. 
This relation may be used to construct a sub-grid model for 
large-scale multi-dimensional simulations, which can in turn 
address how matter and angular momentum are transported 
over the surface. 

In §121 we present the two-dimensional models and demon- 
strate how the mixed mass depends on the velocity of the 
flow. We then, in § |3] incorporate this relation into a one- 
dimensional simulation and discuss two scenarios: that mix- 
ing only occurs if the convective zone reaches the C/O inter- 
face; and that the mixing occurs prior to the onset of con- 
vection. We implement these results as a sub-grid model in 
one-dimensional novae simulations performed with a stellar 
evolution code and discuss the implications, in § |3 for the 
ejected mass of C/O. Because the shear in the envelope prior 
to the onset of convection is unknown, we compute the ac- 
cretion and runaway for several different masses of the mixed 
layers. 

2. BREAKING GRAVITY WAVES 

To calculate the amount of mass mixed by wave breaking as 
a function of time, we performed a suite of numerical simula- 
tions of wind-driven gravity waves. The initial configuration 
consists of two stably stratified layers of compressible fluids 
in hydrostatic equilibrium in a uniform gravitational field g. 
The upper fluid is composed of H/He and the lower fluid is 
composed of C/O. The equation of state is an ideal gamma- 
law of internal energy pe = 1). We choose 7=5/3 
and for the initial condition we construct an isentropic enve- 
lope in hydrostatic equilibrium. The density is then given by 
p(z) = Pi[l - (7- 1)z/hYI^'^^^\ where p,=i is the density 
of the upper fluid just above the interface, p,=2 is the density 
of the lower fluid just below the interface, H ~Y^^cl/g\s, the 
pressure scale height, and is the adiabatic sound speed. This 
choice makes the Brunt-Vaisala frequency zero, i.e., there are 
no internal gravity waves within each of the two layers, which 



isolates the interface modes for study. 
We impose a wind in the upper fluid with velocity 

f/(z)=f/max(l-e"^'/^), (1) 

where z is the vertical direction and z~Q corresponds to the 
location of the initial interface; Umax is the asymptotic maxi- 
mum velocity at z ;2> 5, and 5 is the length scale of the shear 
boundary layer The form of the wind profile (eq. Q) is cho- 
sen to model the boundary layer formed at the surface of the 
white dwarf, and is motivated by studies of the wind-wave 
resonant instability in oceanography (Miles 1957). This in- 
stability is the principal mechanism that drives gravity waves 
and air-water mixing in oceans. In the limit 5 — > 0, one re- 
covers the Kelvin-Helmholtz wind profile. Without a study 
of angular momentum transport in the accreting envelope, we 
do not know the form of the wind profile and in this work 
investi gate the mixing for a range of 5, with 5/// <C 1. 

Alexakis. Young. & RosneJ ^2002^ performed a linear sta- 
bility analysis of the wave amplification for the novae prob- 
lem, and in a forthcoming paper (Alexakis et al., in prepa- 
ration) we will present extensive numerical studies of wave- 
wind interactions. Here we apply those results to the con- 
ditions at the base of the white dwarf envelope. The simula- 
tions a re performed u sing the Flash code (Calder et alJ2000l 
l2?M iFrvxell et al.l ^000 ). a parallel, multi-dimensional, 
adaptive-mesh code for simulating compressible, inviscid 
flow. The computational grid was a uniform mesh of reso- 
lution 1024 X 1024 with periodic boundary conditions along 
the si des of the box and hydrostatic boundaries (Zingal e et alJ 
l2002h at the top and bottom. The interface is perturbed with 
specified modes and a velocity perturbation added to the sur- 
face of the WD representing the motion of the gravity waves. 
The wind profile is specified as an initial condition and in par- 
ticular is not forced. For this study, we set the Mach number 
Ma = t/max/cs — 0.5 and consider four values of 5/H: 0.005, 
0.01, 0.02, and 0.04. These values vary the reciprocal Froude 
number (dg/Ul^^y/^ = {5/H)/{yMa^), which measures the 
available kinetic energy for mixing. 

The unstable modes have wavenumbers jAlexakis et alJ 
"20025 /t>g/f/2^^(l-pi/p2). We chose the size of our com- 
putational domain to be roughly equal to this maximum wave- 
length, so that we have 1024 grid points along the largest pos- 
sible wave. In practice, smaller waves become dominant, with 
the wavelength of the dominant mode depending on the wind 
profile. In all cases this dominant mode is well-resolved. Fig- 
ure [2 shows fully developed waves breaking, generation of 
the "spray", and mixing of the white dwarf substrate up into 
the atmosphere. Mixing is not only effected by the overturn- 
ing of the wave but also by the wind breaking off the tip of 
the wave. Figure|2]shows contours (black lines) of the carbon 
mass fraction at 0.49, 0.20, and 0.02, respectively, for the sim- 
ulation depicted in Figure Q but at a later time. The contour 
at 0.49 corresponds to the carbon mass fraction of the under- 
lying white dwarf, while the smaller contour values indicate 
how far outward into the accreted material the white dwarf 
substrate is mixed. The contour at 0.2 delimits the region that 
contains most of the initial enrichment. 

Figure |3] shows the surface mass density (the mass of C/O 
in the mixed layer per unit area), dMco / dS, averaged over 
the horizontal direction. The mixed layer is defined here as 
the region in which the carbon mass fraction is between 0.49 
and 0.01. The amount of mixed mass depends on these de- 
limiters; reducing the upper limit to 0.4 decreases the mixed 
mass by less than a factor of 2. The C/O is mixed rapidly un- 




Fig. 1. — Breaking CO waves, as determined by simulations in two di- 
mensions. Gravity points towards tlie bottom of the figure, with the vertical 
distance y in units of the pressure scale height H, as evaluated just above the 
interface. The color scale indicates the mass density in units of gem"'. 
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Fig. 2.— Mass fraction of '^C for S/H = 0.04 after t = 3500S/U. The 
vertical dimension is scaled to the pressure scale height H as evaluated just 
above the interface. The contours for '^C mass fractions of, from the top, 
0.02, 0.20, and 0.49. 
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Fig . 3 . — The mixed mass of C/O, per unit area, as a function of time. This 
was computed by averaging over the horizontal direction in the simulations. 
Time is scaled to 5/f/n,ax and Mqo is scaled to piU^^^^/g (see eq. (2)). Four 
different values of gS/U^^ were used, 0.01, 0.02, 0.04, and 0.09. 



Table 1 

Properties of the ID nova models 



Case 


^enrich 




M:nvel 




^nucmax 




(Me) 




(Mo) 




(ergss"') 


no pre-enrichment 







6.93 X 10" 


-5 


6.68 X 10*2 


wave pre-enrichment 


1.26 X 10" 


-<i 


7.08 X 10" 


-i 


7.24 X lO''^ 


wave pre-enrichment 


1.26 X 10" 


-8 


7.05 X 10" 


-5 


7.65 X 10'*^ 


wave pre-enrichment 


1.26 X 10" 


-7 


6.78 X 10" 


-5 


7.71 X 10"^ 


wave pre-emichment 


1.31 X 10" 


-6 


5.42 X 10" 


-5 


1.27 X lO"*' 


wave pre-enrichment 


3.77 X 10" 


-6 


4.67 X 10" 


-5 


2.68 X 10-*' 


wave pre-enrichment 


1.16 X 10" 


-5 


4.85 X 10" 


-5 


1.26 X lO'"* 


enriched accretion 


(50%) 




2.08 X 10" 


-6 


5.28 X lO""' 



til it saturates; further mixing occurs on diffusive timescales. 
Our two-dimensional simulations show, for the range of pa- 
rameters examined, that the total mass of white dwarf mate- 
rial that becomes mixed is independent of the lengthscale 5. 
The rate of mixing, i.e., the initial slopes of the curves in Fig- 
ure |3] does show, however, some dependence on 8. Using 
dimensional analysis and the numerical results, we find that 
for a fixed density ratio pi/p2 = 0.6 (appropriate for the de- 
generate plasma prior to ignition) the total mass per unit area, 
dMco I dS, mixed into the accreted H/He saturates at 



dS 2 



(2) 



where a is a non-dimensional constant' ' that we determine 
from the simulations to be a ~ 0.6 (see Figure |3}- The 
timescale to reach saturation is far shorter than the timescale 
of either the accretion phase (> 10^ yr) or the pre-peak con- 
vective phase (~ 100 yr) of a typical classical nova. For exam- 
ple, if 5 /H ~ 0.01 and Ma = 0.5, then the saturation timescale 
is of order seconds for our nova setup. 

3. one-dimensional NOVA MODELS 

We now incorporate our simulations of the wave break- 
ing and mixing into simulations of the thermonuclear 

" In general, a depends on P2/P1 and possibly Ma. A parameter study, 
outside the scope of this paper, is necessary to determine this dependence. 



burning of a nova. In order to explore the global 
properties of this local mixing mechanism, we com- 
pute several one-dimensional models of novae with a 
modifi ed version of t he KEPLER stellar evolution code 
JWeaver. Zimmerman. & Wooslev 1978) that includes a large 
network of ~ 100 light isotopes (R auscher et al.^2002.) implic- 
itly coupled into the stellar structure solver. This allows us to 
follow throughout the entire convective envelope the radioac- 
tive decays of nuclei that are formed at the base of the burning 
zone. The underlying white dwarf has a mass of I.OMq, a 
radius of 5000km, a luminosity of 10^'ergss^', and is com- 
posed of a 50%/50% '^C/'^O mixture. We retain the outer 
O.OO5M0 on the computational grid to follow the thermal in- 
ertia of the outer WD layers and allow the model to relax un- 
til it has a constant luminosity before we start the accretion. 
The accretion rate is lO^^Mp^yr^', and in all but one model 
(see item 3 below), the accreted material consists of 70 % 'H, 
28.7 % ^He, 0.3 % '^c, 0.1 % ''^N, and 0.9 % '^O, by mass. 
Convection is modeled using the Ledoux criterion for stabil- 
ity and mixing length theory. We assume no convective over- 
shooting and that semiconvection is too slow to cause mixing. 

We investigate three scenarios, as summarized in Table ^ 
for generating the gravity wave induced mixed layer 

1 . In the first scenario, the only shearing considered is that 
from the convective cells driving a wind at the inter- 
face between the H-rich atmosphere and the C/O sub- 
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Fig. 4. — Kippenhahn diagram of a nova without enrichment. The ,t-axis 
indicates time intervals for the different evolution stages, and the y-axis gives 
the mass above the C/O WD substrate. Green hatching (framed by a green 
line) indicates convection, blue shading indicates nuclear energy generation 
for which each level of darker blue denotes an increase by one order of mag- 
nitude, starting at 10'" ergs g^' s^' . The thick black line shows the total mass 
of star (including ejecta), increasing because of accretion; the dash-dotted 
line indicates the mass outside of 10'~cm; and the dashed line marks the 
interface between the white dwarf C/O substrate and the accreted layers. 



strate. This is the scenario envisaged by ' Rosner et alJ 
(2001), in which the C/O, after being entrained, is then 
distributed throughout the convective zone. Figure 0] 
summarizes the evolution of the accreted layer. After 
accretion for about 6.9 x 10^ yr a convective zone forms 
about 2.5 X lO^^M© above the WD interface. The total 
envelope mass accreted at ignition is w 7 x 10^^ 
in good agreement with the estimates of Fuiim otol 
dllllS. Within - lOOyr from the onset of convec- 
tion, the convective zone extends upward to the surface 
and downward to about 5 x 10^*' above the inter- 
face, at which time the runaway reaches its peak rate 
of energy generation. This peak phase evolves on a 
timescale of ~ 1 h, but the convection does not reach 
the WD interface. Only days later, when the nova enve- 
lope is already significantly expanding, does the burn- 
ing layer reach the WD. This downward movement, in 
mass, of the burning layer is mediated by heat conduc- 
tion into these deeper H-rich layers. The convective 
zone also moves deeper, but never quite reaches the WD 
interface^^. As a result, there is no injection of C/O into 
the H-rich envelope, by construction. 

The envelope expands to large radii (> 10^^ cm, large 
enough to engulf the secondary) by the end of the sim- 
ulation, but the result will only be a slow nova, in 
ro ugh agreement with the semi-analytical calculations 
of lRiiimotol Jl982albft . We do not find any contact of 
the convective region with the WD substrate; as a re- 
sult, we do not expect that any convection-driven wave 
mixing would occur, unless there is a large redistri- 
bution of angular momentum in the accreted envelope 

At the very bottom of the H shell the nuclear luminosity from H burn- 
ing can never become large enough to drive convection, as it goes to zero at 
the interface; the only possibility for mixing to occur is if a steep tempera- 
ture gradient arises from the thermal inertia of the heated WD core below an 
expanding H envelope. 




C/O core - 



Fig. 5. — Same as Fig.l^but for a model in which the inner M„,ix = 4.6 x 
10^* g are enriched in C/O with a linear composition gradient (with respect 
to the Lagrangian mass coordinate) between the WD composition (C/O) at 
the base and the accretion composition (solar) at the upper edge. Note that 
the convective zone does not reach the interface with the WD substrate, and 
that a significant semi-convective region, indicated by red hatching, develops 
prior to the onset of convection. 



(iKutter & Sparksll989tlSparks & Kutteill987h . 

. In the second scenario, the shearing originates from 
a wind at the base of the accreted envelope blowing 
across the underlying C/O substrate. We assume that 
the wind persists throughout the H/He layer with veloc- 
ity sufficient to drive mixing on a timescale much less 
than that to accrete a critical mass of fuel. In this case 
the mixed layer is generated prior to ignition. We also 
assume a linear (in Lagrangian mass coordinate) gradi- 
ent in the mass fraction of C/O between the WD and the 
accreted envelope. Because the shear profile, and hence 
the amount of mass mixed, is unknown, we consider a 
range of mixed masses (see TableQ. Our ignorance of 
the fluid motions in the envelope and substrate prevents 
us from saying where and how the accreted material 
spreads over the surface, and so we cannot determine 
the actual shear at the base of the envelope. In all cases 
that we consider, the shear velocity is subsonic, which 
is the range of validity of the 2-d simulations. For such 
velocities, the mixed layer is always thinner than a pres- 
sure scale height (eq. ||2l)- 

Figure |5] shows the case with the largest pre- 
enrichment: the mixed material comprises 4.6 x 10^^ g 
(i.e., a total of 2.3 x 10^^ g of WD material is being 
mixed with H-rich material). This is about 25% of 
the envelope mass, which is less than that generated 
by the calculation described in section |2] From equa- 
tion (|2j, a 25% enrichment corresponds to Ma = 0.4 
if the velocity were uniform over the surface. We 
note that this velocity is much less than Keplerian, 
f^max ~ 0.05(GM//?)'/^, as expected in the envelope 
well below the accretion disk boundary layer (see, e.g., 
Ifopham & Naravan 1995). If the velocity were not uni- 
form over the surface, then a higher maximum velocity 
would be required to inject the same percentage of the 
total envelope mass. 

The mixed layers are added to the surface using 
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Fig. 6. — Time-integtated average specific nuclear energy generation of 
the accreted enevelope. Tlie curves are shifted so that the zero in time cor- 
responds to the peak energy generation rate. The gray curve indicates the 
model with the highest wave-induced pre-enrichment (see Tab.fTland Fig. 151; 
the black curve indicates the model without prior enrichemnt (Fig.|3. The 
slope of a curve is proportional to the average specific energy generation rate 
in the envelope. 



the same accretion rate as the rest of the envelope 
(lO^^M0yr^'); to prevent a runaway while these lay- 
ers are being added, we suspend nuclear energy gener- 
ation during the accretion of the first 5 x 10^** g. About 
50 yr prior to peak energy generation a semiconvective 
region forms in the outer 60 % (by mass) of the enriched 
layer; about 8 months before the peak a convection 
zone starts at the base of this semiconvective layer (at 
7 X 10"*' Mq above the WD interface) that extends up- 
ward to the surface and downward to about 5 x IO^^Mq 
above the interface. In contrast to the first case, the pre- 
mixed envelope with the largest enriched mass ignites 
at a smaller accumulated mass, « 5 x 10-^ Mr.,. 

The energy generation in the enriched layer is domi- 
nated by i2c(p,7)i3N(j3 + )i3c. The peak' energy gen- 
eration timescale is about 3 min — 20 times faster than 
in the first case — with a specific energy generation rate 
about 20 times higher. Figure |6l shows the specific nu- 
clear energy generation in the accreted envelope, with 
the time centered about the peak energy generation 
rate. The larger rate of energy generation roughly cor- 
responds to the ^ 20 times greater metal enrichment. In 
this scenario the entire 25 % of pre-enriched C/0 mate- 
rial is spread throughout the hydrogen envelope during 
the runaway. This leads to a a corresponding enrich- 
ment in the ejecta. The ejecta are unbound, i.e., they 
have a positive velocity at very large radii. 

As in the first case, the convection never reaches the 
WD interface (although the base of the convective zone 
is much deeper), so there is no additional convection- 
induced wave mixing of the C/O substrate with the 
hydrogen-rich envelope. We did experiment with pre- 
runaway, pre-enriched layers as small as 1/10,000 the 
mass used in the case described above (TableQ. As be- 
fore, we disregarded nuclear burning for the first 10^^ g 
of accreted material, with the exception of the case with 
3.77 X 10^^ Mq of enrichment, where we took twice 



that mass to accommodate the entire enriched layer. 
The difference in masses for which nuclear burning is 
suspended is the reason for the case with the largest pre- 
enriched mass, 1.15 X lO^^Mp having a slightly larger 
envelope mass at runaway than that of the case with 
3.77 X lO^^M© of enriched mass. 

In none of these cases did the convection reach the 
WD interface. Eventually, for the case of an enriched 
mass of < lO^^M©, the runaway behaves as in the non- 
enriched case: the energy release rate in the mixed lay- 
ers is too small to drive convection above and heat con- 
duction cools it efficiently so that the runaway does 
not occur in its vicinity. The slightly higher enve- 
lope masses obtained in the limit of small enrichment 
is likely to be an artifact of our suspension of nuclear 
energy generation for the first fraction of the accretion 
phase. 

3. In a third scenario we assume accretion of material en- 
riched to a mass fraction of 50 % in C/O. Allowing for 
nuclear energy generation instantaneously in this mate- 
rial leads to a rapid runaway after only a thin layer has 
been accumulated at the surface of the star. If this mag- 
nitude of enrichment was due to instantaneous mixing 
with a colder WD substrate, then a lower temperature in 
the enriched material and thus a significantly later run- 
away could result. We include this run to demonstrate 
how qualitatively different the accretion of pre-enriched 
material behaves as compared to first accreting unen- 
riched material and then, only later after it has settled, 
mixing the accreted layers with the C/O substrate. 

4. IMPLICATIONS AND SUMMARY 

Using a constant wind profile that blows across the sur- 
face of a white dwarf, we have performed a two-dimensional 
param eter study of the mechanism proposed by Rosner et al] 
J2001i) . Our primary results of the mixing rate and the max- 
imum mixed mass, equation (|2ji, suggest that this process 
can mix about lO^^-lO^^M© of the underlying C/O into the 
hydrogen-rich envelope. From this, we investigated two sce- 
narios: that the wind is the result of convection, and that the 
wind originates during the accretion phase. 

We find that if no enrichment occurs prior to the onset 
of convection, then the convective zone does not reach the 
C/O interface, and no additional mixing occurs (in the one- 
dimensional model) in the absence of convective overshoot. 
The result in this case will be a slow nova, with little enrich- 
ment of the ejecta. In contrast, an envelope with a mixed layer 
at the C/O interface, consistent with the scalings from high- 
resolution numerical studies, provides a more violent runaway 
and the ejecta are enriched by 25% in C/O, consistent the 
enrichment observed in some nova ejecta. Such an event does 
require a strong shear velocity, but the saturation amount of 
mass mixed is roughly independent of the shape of the shear 
profile, so long as its thickness is much less than a pressure 
scale height. Our runs, outlined in Table ^ show that when 
the metals are concentrated at the base of the accreted enve- 
lope, the amount of mass accreted prior to runaway is larger 
than if the C/O were uniformly distributed over the envelope, 
i.e., if the white dwarf were to accrete material with a superso- 
lar metallicity. The reason is that the opacity of the envelope 
is larger than when the C/O are concentrated at the base. 

There are a number of issues that we have not yet addressed 
or investigated. First, our results might depend on the dimen- 
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sionality of the system. We would therefore expect that al- 
though the generation of the gravity waves is captured by two- 
dimensional dynamics, the energy cascade of the waves and 
the advection of spray and vorticity will be different in a three 
dimensions. Second, this work investigates only the density 
ratio for an accreted envelope in thermal equilibrium. During 
the early phase of accretion and during the runaway the den- 
sity ratio will likely be different. Third, although our param- 
eter study covers more than one order of magnitude in 8, it is 
important to know if our results still hold at even larger val- 
ues. As 5 becomes a respectable fraction of H, we expect that 
the stratification of the envelope, and in particular the effect of 
a non-zero Brunt-Vaisala frequency, will become important. 
All of these differences can affect the mass and thickness of 
the mixed layer, and the investigation of these issues is the 
subject of ongoing work. 

In our one-dimensional nova simulations, we assumed that 
a convective roll will only entrain the mixed layer if the base 
of the convective zone (as computed from a mixing-length 
formalism) reaches the interface. Future multidimensional 
studies can inform us as to how convection interacts with the 
interfacial gravity waves, and to whether this assumption is 
in fact correct. Moreover, our sub-grid model, when incor- 
porated into a one-dimensional calculation, does not account 



for spatial variations of the wind. Such a variation would oc- 
cur, for example, if the fuel accretes non-uniformly over the 
surface of the star and then spreads. In reality the amount of 
mass mixed will vary over the surface of the star Our one- 
dimensional calculations must of necessity take the mixed 
mass as a free parameter. Because the amount of C/O en- 
trained depends quadratically on the velocity in the wind pro- 
file, our calculation underscores the need for detailed simula- 
tions of the shear profile of an accreting white dwarf These 
issues are clearly important and are the next steps to pursue. 



We thank Lars Bildsten and Stan Woosley for helpful dis- 
cussions. This work is supported in part by the U.S. De- 
partment of Energy under Grant No. B341495 to the Center 
for Astrophysical Thermonuclear Flashes at the University of 
Chicago. A. H. is supported by the Department of Energy un- 
der contract W-7405-ENG-36. L. J. D. acknowledges support 
by the Krell Institute CSGF M. Z. is supported by DOE grant 
number DE-FC02-01ER41 176 to the Supernova Science Cen- 
ter/UCSC. P. M. R. acknowledges support from the Univer- 
sity of Illinois at Urbana-Champaign and the National Cen- 
ter for Supercomputing Applications (NCSA). K. Olson ac- 
knowledges partial support from NASA grant NAS5-28524. 



REFERENCES 



Alexakis, A., Young, Y.-N., & Rosner, R. 2002, Phys. Rev. E, 65, 026313 
Calder, A. C, et al. 2000, in Proc. Supercom puting j OOO, (IEEE Computer 

Soc.) http://sc2000. org/procee dings /| 
Calder, A. C, et al. 2002, ApJS, 143, 201 
Fryxell, B., et al. 2000, ApJS, 131, 273 
Fujimoto, M. Y. 1982a, ApJ, 257, 752 
Fujimoto, M. Y. 1982b. ApJ, 257, 767 
Fujimoto, M.Y 1988, A&A, 198, 163 

Glasner, S.A., Livne, E., & Truran. J.W. 1997, ApJ, 475, 754 
Gehrz, R. D., Truran, J.W., Williams, R.E., & StaiTlield, S. 1998, PASP, 110, 
3 

Hemanz, M., Jose, J., Coc, A., & Isern, J. 1996, ApJ, 465, L27 
Kercek, A., Hillebrandt, W., & Truran, J.W. 1998, A&A, 337, 370 
Kercek, A., Hillebrandt, W., & Truran, J.W. 1999, A&A, 345, 831 
Kippenhahn, R. & Thomas, H.-C. 1978, A&A, 63, 625 
Kovetz, A. & Prialnik, D. 1985, ApJ, 291, 812 
Kutter, G.S. & Sparks, W.M. 1989, ApJ. 340. 985 

Livio, M., & Truran, J. W. 1990, in Nonlinear Astrophysical Fluid Dynamics, 
ed. J. R. Buchler & S. T. Gottesman, (New York: New York Academy of 
Sciences), 126 

Livio, M. & Truran, J.W. 1994, ApJ, 425, 797 



MacDonald, J. 1983, ApJ, 273, 289 
Miles, J. 1957, J. Fluid Mech., 3, 185 
Popham, R. & Narayan, R. 1995, ApJ, 442, 337 
Prialnik, D. & Kovetz, A. 1984, ApJ, 281, 367 

Rauscher, T., Heger, A., Hoffman, R. D., Woosley, S. E. 2002, ApJ, 576, 323 
Rosner, R., Alexakis, A., Young, Y.-N., Truran, J. W., & Hillebrandt, W. 2001, 

ApJ, 562, LI 77 
Sparks, W. M., & Kutter, G. S. 1987, ApJ, 321, 394 

Starrfield, S.G., Truran, J.W., Sparks, W.M., & Kutter, G.S. 1972, ApJ, 176, 
169 

Starrfield, S.G., Truran, J.W., Wiescher. M.C., & Sparks. W.M. 1998, 
MNRAS, 296, 502 

Truran, J. W. 1982, in Essays in Nuclear Astrophysics, ed. C. A. Barnes, D. D. 

Clayton, & D. N. Schramm (Cambridge: Cambridge University Press), 467 
Weaver, T. A., Zimmerman, B.A., & Woosley, S.E. 1978, ApJ, 225, 1021 
Woosley, S.E. 1986, in Nucleosynthesis and Chemical Evolution, Eds. 

B. Hauck, A. Maeder, & G. Magnet (Savverng, Switzerland: Geneva 

Observatory Publication), p. 1 
Zingale, M., et al. 2002, ApJS, 143, 539 



